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Coarse-graining atomic displacements in a solid produces both local affine strains and “non-affine” 
fluctuations. Here we study the equilibrium dynamics of these coarse grained quantities to obtain 
space-time dependent correlation functions. We show how a subset of these thermally excited, 
non-affine fluctuations act as precursors for the nucleation of lattice defects and suggest how defect 
probabilities may be altered by an experimentally realisable “external” field conjugate to the global 
non-affinity parameter. Our results are amenable to verification in experiments on colloidal crystals 
using commonly available holographic laser tweezer and video microscopy techniques, and may lead 
to simple ways of controlling the defect density of a colloidal solid. 


I. INTRODUCTION 

While a large body of work has accumulated over many 
decades on the physics of crystal defects [T], the micro¬ 
scopic causes of defect nucleation and yielding in solids 
remain active areas of recent research mm- A small 
external stress on a crystalline solid at non-zero temper¬ 
atures affects atomic configurations in two ways: (1) an 
affine deformation characterised by the elastic strain and 
(2) a modification of the relative probabilities of ther¬ 
mally excited lattice defects ms]. Within a linear re¬ 
sponse picture [5] local strain fluctuations measured from 
particle coordinates determine the elastic moduli of the 
solid, which in turn govern the magnitude of the affine 
response EHin]. It is therefore natural to ask the com¬ 
plementary question viz. fluctuations of which quantity, 
derivable solely from the configuration of the atoms, mea¬ 
sure the susceptibility of a crystalline solid to creation of 
defects? In this paper, we pursue this issue by extending 
and generalising an approach introduced in m based on 
coarse-graining of atomic displacements. Soft, precursor 
fluctuations which give rise to defects appear as a natu¬ 
ral outcome of this coarse-graining process. We explore 
some of the interesting consequences of this connection 
- such as the ability to engineer equilibrium defect con¬ 
centrations, at least in a colloidal crystal na, by subtly 
altering the statistical weights of these precursors using 
dynamic light fields [13]. 

Consider a system consisting of i = 1.. .N particles 
with instantaneous positions {r} vibrating about a set 
of reference coordinates {R}. To begin, we first elevate 
a measure of non-affinity introduced in El to identify 
elastic heterogeneities in sheared amorphous solids, to 
the role of a fully-fledged thermodynamic, collective co¬ 
ordinate. This variable, X = N~^ ^ scalar 


functional of both the instantaneous and the reference 
coordinates, measures the magnitude of non-elastic devi¬ 
ations of the positions of all particles away from their 
reference configuration coarse-grained over a reference 
volume U. The local x(Ri) is a function of the instanta¬ 
neous and reference positions r and R of the particles in 
the neighborhood U of a given particle i with reference 
particle position R^. We had earlier obtained the equilib¬ 
rium statistics of of y (spatial dependence suppressed for 
brevity), in crystals [11] at finite temperatures. We had 
shown that under an external stress 5], particles undergo 
both affine and non-affine deviations, with y always in¬ 
creasing as within the harmonic approximation. 

Here, we first extend the work reported in El to time- 
dependent correlation functions for y and strains at zero 
stress. This part of our treatment is similar in spirit 
to that of [15] where the dynamical density correlations 
are analysed in terms of a sum over harmonic degrees 
of freedom. The relaxation time of an observable arises, 
within this approximation, from the de-phasing of in¬ 
coherent harmonic oscillations. While our treatment is 
perfectly general and is applicable to any solid in any di¬ 
mension for which {R} and the interactions are known, 
we present results for the two dimensional triangular lat¬ 
tice. Next, an analysis of the vibrational modes con¬ 
tributing to non-affine distortions of U reveals that most 
of X arises from two degenerate non-affine displacements 
which tend to replace four 6 coordinated particles with 
two pairs of particles with 5 and 7 neighbours: an incipi¬ 
ent, or precursor, dislocation-anti-dislocation pair. Hav¬ 
ing obtained time dependent correlation functions of x, 
we identify a field, hx^ thermodynamically conjugate to 
the global non-affinity parameter X and argue that the 
response of the solid to an application of this field would 
change the magnitude of X (and x) and hence the prob¬ 
ability of these localised precursor fluctuations that give 
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rise to dislocation pairs. Since hx couples linearly, X can 
both increase or decrease depending on the sign of hx- 
Further, unlike stress, hx changes defect probabilities 
directly without introducing affine strains and associated 
spatial anisotropies, simplifying the study of their statics 
and dynamics. Finally, since X is given entirely in terms 
of the particle coordinates and the {R} which appear as 
constant parameters, hx can be directly introduced into 
the Hamiltonian and the dynamics of system calculated 
using standard molecular dynamics m- For a colloidal 
solid, it is even possible to apply hx in the laboratory us¬ 
ing laser tweezers [HI HU, a point that we elaborate later. 
We explore the effect of hx on the equilibrium statistics 
and dynamics of X(and y) and of the individual defect 
precursor modes. 

Before we end this introduction, we mention an aspect 
of our work that makes it relevant to some recent studies 
of the mechanical response of soft amorphous solids [18] 
and glasses [19]. In such solids, it is impossible to define 
the kinds of defect configurations encountered in crystals 
such as vacancies, dislocations, stacking faults or grain 
boundaries [20]. However, given any reference configu¬ 
ration and a set of particle coordinates, y and much of 
everything else discussed in our work can still be defined 
and computed. The precursor modes in this case should 
be related to non-affine droplet fluctuations, which have 
been extensively studied in recent years [HllIHlg]. Un¬ 
like crystals, however, the identification of defect precur¬ 
sors with actual defects is much more problematic for 
amorphous solids due to the lack of a simple and unique 
reference {R} and is, therefore, a subject of ongoing 
lively debate [30H32]. We hope that some of the ideas 
discussed here may be useful in illuminating this issue. 
We return to a discussion of this point briefly later. 

The rest of the paper is organised as follows. In sec¬ 
tion [H] we set up the calculation and define the coarse- 
graining process used to calculate spatio-temporal corre¬ 
lation functions. Parts of this calculation have previously 
appeared in m , but we include the relevant aspects here 
for completeness and to make the paper self-contained. 
In [lll|we present our results for the time dependent y 
and strain fluctuations in the two-dimensional triangular 
lattice. In [^we identify defect precursors and obtain 
their statistics. We also introduce the non-affine field hx 
and study its effect on these precursor fluctuations. In 
sectionjVjwe suggest how hx may be produced in the lab¬ 
oratory using laser tweezers m- We discuss our results 
and conclude by giving indications of future directions in 
section ED 


II. COARSE GRAINING AND DYNAMIC 
CORRELATION FUNCTIONS 

On application of an external stress or as a result of 
thermal fluctuations, particles i within a solid undergo 
displacements — R^ away from some chosen ref¬ 

erence configuration R^ to their displaced positions r^. 


In a homogeneous solid at vanishing temperature, such 
displacements are affine, implying that they can be ex¬ 
pressed as = DR^, where D = is the defor¬ 

mation tensor related to the external stress 5] via the 
tensor of elastic constants K. To derive the closest ap¬ 
proximation to this simple zero temperature scenario in 
the presence of thermal fluctuations we proceed as follows 

m- 

Consider a neighbourhood, O, larger than the unit cell, 
around a central particle labelled 0 consisting of Nq par¬ 
ticles i within a cut-off distance Rq in a d dimensional 
lattice. The reference, zero temperature lattice config¬ 
urations are labelled by Ri=o...Arj^ while the fluctuating 
atom positions are u=o...Ar5nj- The particle displacements 
are then as before — R^. Now define relative 

displacements, — uq = — tq — (Ri — Ro) of 

particle i compared to particle 0. The “best fit” [14] 
coarse-grained local deformation tensor D is the one that 
minimizes ~ D(R^ — Ro)]^ with the non-affinity 

parameter y being the (positive definite) minimum value 
of this quantity. 

In m we showed that the result of this minimi¬ 
sation procedure may be expressed as a projection of 
the particle displacements into mutually orthogo¬ 
nal subspaces as defined by two projection operators 
P and RQ. In terms of these, y = A^PA while the 
elements of the affine deformation tensor (strains and 
local rotation), Dq,^, arranged as a linear array e = 
(Dii,Di 2 , ... ,Di^,D 2 i, .. .,Ddd), are given by e = QA. 
Here A is a column vector with Nd elements contain¬ 
ing the components of the A^. The projectors are given 
explicitly by RQ = R(R^R)“^R^ and P = I — RQ. As 
in m, the Nd X d? matrix R appearing here has ele¬ 
ments Ria, 77 ' = ^aj{Ri^' — Rqy) where the Riy and 
Ro, 7 ' are the components of the lattice positions R^ and 
Ro, respectively. Now define the correlation matrix C 
with elements, Ciaj^y = {Ai^Aj^) where the angular 
brackets (...) indicate an average over the equilibrium 
ensemble. One can then easily obtain the statistics of y 
and e in terms of C. For example the probability dis¬ 
tribution for the affine distortions e is a dimensional 
Gaussian with zero mean and co-variance matrix QCQ^ 
whose elements are proportional to the elastic moduli. 
On the other hand, y is distributed as the sum of the 
squares of NQd — d? independent Gaussian random vari¬ 
ables with variances given by the eigenvalues of PCP. A 
comparison of the projected atomic displacements, i.e. 
eigenvectors of PCP and (I — P)C(I — P), that give rise 
to the y and e shows that while the latter consist of lo¬ 
cal volume, uniaxial and shear distortions of Q together 
with local rotations, non-affine displacements, which con¬ 
tribute to y, correspond to small wavelength distortions 
of particles within Q. Application of an external stress, 
X), shifts the strain probability distributions to non-zero 
mean strain in accordance with Hooke’s law and fluctu¬ 
ation response relations but does not affect y to linear 
order. The lowest order variation of y with 5] is given 
by (x)s = (x)s=o + S'^QC[P, C]Q'^S, where [P, C] is a 
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commutator. 

In order to calculate the spatio-temporal correlation 
functions of the non-affinity y and strains e, we need to 
consider simultaneously displacement differences in two 
neighborhoods O and Cl centered on lattice positions Rq 
and Ro at time t and t' respectively. The vector A(t) is 
defined as the displacement corresponding to the refer¬ 
ence lattice position Rq at time t, with an analogous defi¬ 
nition for A(t'). The local afhne strain e(Ro, t) = QA(t) 
and non-affinity x(Ro,t) = A^(t)PA(t) are defined as 
before. For time t' and position Rq we have the corre¬ 
sponding quantities e(Ro,t') = QA(t') and x(Ro,t') = 
A^(t')PA(t'). The covariances may now be defined as 

CiaJj = = (A^Q,(0)Aj^(0)), 

CiaJ^ = (A(t )ia^{t )jj) = (A(O)ic^A(O)j^) 

( 1 ) 


interactions we may approximate as being harmonic. Al¬ 
ternately, the covariance matrix may also be obtained 
experimentally [9] in the case of colloidal solids using 
video microscopy without any a priori assumption con¬ 
cerning the form of the interactions. One may directly 
measure (uqU^q,) = 'D~^ {q) vbz^{^. — ^ where T>(q) is 

the dynamical matrix, and vbz the volume of the Bril- 
louin zone. Given the dynamical matrix, Cia^j'y may 
be evaluated as follows. We substitute for the rela¬ 
tive displacements their expansion in terms of the vi¬ 
brational modes of the lattice viz., A^(t) = — uq = 

^'^BzEs/c^quTas(q)a5(q)(e*‘»-*^» - cos[w5(q)t], 

into the third of the equations Q to obtain, 

/=v _ /" j ^ / .cos[ws(q)t] 

Cia,j^ — /. I ( q ) ( q ) 2^ \ ^ 

vbz^J 

(4) 


Obviously the first two averages are identical and reduce 
to the space and time independent second-order moments 
(AAT) [TT]; the third quantity yields the required corre¬ 
lation functions. To derive the expressions for the time- 
dependent strain and non-afhnity auto-correlation func¬ 
tions we use their definitions in terms of the relative dis¬ 
placement projections. We obtain, therefore, 

Ce(Ro,i,Ro7') = (e(Ro,f)e'’’(Ro,f')) 

= (QA(f)A'^(f')Q^) 

= QCQT. (2) 

The correlation functions between any pair of affine 
strains may now be obtained by taking appropriate linear 
combinations of the elements of Ce- In the next section 
we focus on one such component, viz, the shear strain e. 
Similarly, the correlation between x(Ro 5 ^) and 
can be calculated using Wick’s theorem as 

C'x(R-o,f,Ro,i') = (x(Ro7)x(Ro,f')) - ixf 

= 2Tr(PCP)(PCP)^ = 2^d| (3) 
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where, in the final equation, the denote the Nq^ d — dP 
non-zero eigenvalues of the matrix (PCP)(PCP)^. Of 
course, in a homogeneous solid in equilibrium, these cor¬ 
relation functions are functions only of the relative coor¬ 
dinates Ro — Ro and times t — t'. We will denote these 
simply by R and t in what follows. 

Note that so far we have not made any assumptions 
about the structure and interactions of the particles i 
and all our results apply equally well for any system in 
any dimension as long as a well defined reference config¬ 
uration {R} exists. Indeed, we believe that a fair frac¬ 
tion of our results should apply even to amorphous solids 
with displacements being measured from a set of par¬ 
ticle coordinates obtained from a zero temperature en¬ 
ergy minimisation. To obtain analytic results we need 
to evaluate the covariances and for the rest of this pa¬ 
per we specialise to periodic lattices of particles, whose 


In the above expressions I is the lattice parameter and 
as(q) and 0Js{^) ^.re the eigenvectors and eigenvalues 
(phonon frequencies) respectively of the dynamical ma¬ 
trix corresponding to the phonon branch. The q- 
space integrals are over the Brillouin zone. 


III. RESULTS FOR THE 2D TRIANGULAR 
CRYSTAL 


The formulation for the spatio-temporal correlation 
functions given in the previous section (section is ap¬ 
plicable for any periodic crystal as long as the dynamical 
matrix is known. In this section we present our re¬ 
sults for the simple but important case of a triangular 
network of particles connected by harmonic springs de¬ 
fined by the Hamiltonian, 




harm 



(b‘) 


( 5 ) 


where u^, and m are displacement, momentum and 
mass of the particle i respectively, and the unit vector R^j 
points from particle i to j in the reference configuration. 
The sum in the second term in ^ runs over all bonds in 
the network, each with spring constant K. The unit of 
distance will be the lattice parameter from now on while 
time will be measured in units of K. The tempera¬ 
ture may also be rescaled to unity without loss of general¬ 
ity. Because of its simplicity, the harmonic triangular net 
has been studied extensively and is known to be a good 
approximation for many real crystalline solids in two di¬ 
mensions HHIO]. The dynamical matrix and hence the 
dispersion relation co’(q) for this system is also known [33]; 
for small q it is given by 0Js=t,l = Cs=t,l\^\ with the 

transverse and longitudinal sound velocities ct = ^ 

and Cl = |y^- We consider a coarse-graining volume 
Cl consisting of a central atom and its Nq = 6 nearest 
neighbours in the triangular lattice m- The normalised 
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FIG. 1. (a) Cx(0,t) from computation (line) compared with 
that obtained from molecular dynamics simulations (points) 
of a 500 X 500 site harmonic lattice with unit particle masses. 
Note that the error bars are smaller than the size of the sym¬ 
bols. (b) Plot of Ce(0,t) for the same system as in (a). 


correlation functions for x and e fluctuations are given 
by 


and 




(x(o,o)x(R-.^)) - (x)^ 

(x^) - (x)^ 


Ce{R,t) 


(e(0,0)e(R,i)) 

(e2) 


where e = (1^12 + D 2 i )/2 = es is the shear strain. Note 
that the temperature (or coupling constant) scales out 
for these quantities. 

In Figg a), we plot C^(0,t). The integrals over the 
Brillouin zone were computed numerically using a 256 
point Gauss quadrature for each real space lattice posi¬ 
tion. Our results are compared with those obtained from 
molecular dynamics simulations of a 500 x 500 site har¬ 
monic triangular net. To obtain the simulation results, 
the solid was first allowed to equilibrate for 4 x 10^ MD 
steps with a time-step of 5 x 10“^ in the canonical (NVE) 
ensemble. Data was collected from 4 x 10^ to 4.2 x 10^ 
MD steps at intervals of 400 MD steps. Correlation func¬ 
tions were obtained by averaging over particles. 

It is clear that our results agree with simulation data 
within the error bars of the latter. The decay of is 
not monotonic but oscillatory, a feature arising from the 
time-periodic lattice vibrations of the solid, which are 
all in phase at t = 0, gradually de-cohering for larger 


times m- Similar oscillations are also observed in 
shown in Figj^b) for the same system and both x and e 
relax over similar time scales. 

The complex relaxation of the dynamic correlations 
is even more in evidence when we evaluate these func¬ 
tions in both space and time. In Figj^a)-(c) we plot 
the full (R, t) for the first few nearest neighbour lat¬ 
tice points of the 2d triangular lattice. For small times, 
this correlation function is sharply peaked at the origin 
and decays rapidly to zero after the second neighbour 
shell (FigJ^a)). At larger times, the function decays but 
becomes longer ranged extending up to the sixth neigh¬ 
bour shell for t = 5 (FigJ^c)), corresponding to a spread 
with a speed comparable to ct] finally, C^(R,G ^ 0 
everywhere for large time differences t. In Fig. [^a)-(c) 



FIG. 2. Cx(R, t) calculated at different times (a) t = 0, (b) 
t = 2, (c) t = 5 for the harmonic triangular lattice. The 
parameters are the same as in Fig[^ Note changes in scale 
between the three subpanels. 


we show the corresponding space-time correlation func¬ 
tions for the shear strain. The equal time spatial corre¬ 
lation FigJ^a) has been calculated before m and has 
also been measured from video microscopy of colloidal 
solids [To]. The typical four-fold symmetry (butterfly pat¬ 
tern) of this correlation function has also been observed 
in experiments on amorphous colloids m- This pattern 
is easy to understand since (7e(R, 0) represents the re¬ 
sponse of a solid to a delta function shear load at the ori¬ 
gin, which may arise from a small inclusion or “Eshelby” 
defect [22]. At subsequent times, the correlation func¬ 
tion retains its significant four-fold symmetry, although 
it shows wave-like oscillations in space and time; these 
eventually decay to zero. 

The four-fold symmetry of the strain correlation func¬ 
tion can also be understood from a momentum space 
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FIG. 3. Ce(Il,t) calculated for times (a) t — 0, (b) t = 2, 
(c) t = 5 for the same set of parameters as in Fig.m Note 
changes in scale between subpanels. 


formulation as described in m , where the Fourier trans¬ 
form of the shear strain correlation function is simply 
given by 

= Ei 212 + E212I + 2E122I, 


with 

Eaa'77'(q, t) = ^ (g) 


Substituting the expressions for E into and expand¬ 
ing all q dependent quantities to leading order in the 
wavenumber, we finally get (g = |q|) 


Ceiq,t) 


‘^QxQv 

2 J COs(CLg^) + 


c|,q4 


cos{ct qt). 


The four-fold symmetry of is now obvious. 


IV. DEFECT PRECURSORS IN THE 2D 
TRIANGULAR CRYSTAL 

In the previous two sections we derived a systematic 
procedure for analysing particle displacements within a 
coarse-graining volume Q as affine or non-affine. The 
affine displacements can be identified as elastic strains 
whose fluctuations determine the elastic constants of the 
solid. In this section we turn to the identity of the non- 
affine fluctuations and show that the fluctuations with 
the highest contribution to y represent precursors to 
the formation of pairs of lattice defects. We also show 
that they are statistical fluctuations that obey standard 


fluctuation-response relations and thereby identify the 
conjugate field hx- Positive values of hx enhance and 
negative values suppress lattice defects. Finally we cal¬ 
culate space-time correlation functions for these lattice 
distortions in the presence of nonzero hx^ using results 
derived in the earlier sections. 

Recall that in the 2d triangular lattice (y) = 
where the are the eight non-zero eigenvalues of the 
Nq d X Nq d = 12 X 12-dimensional matrix PCP HD. 
The eigenvectors corresponding to these eigenvalues 
represent non-affine distortions of the coarse graining vol¬ 
ume, their relative contributions to y being determined 
by the value of In Fig. we plot the magnitudes of 
It is immediately clear that there are three groups of 
terms. The eigenvalues of the two degenerate, non-affine 
modes corresponding to /i = 1 and ji = 2 are separated 
from the next higher one /i = 3 by a large gap - a fac¬ 
tor of 4 - and from the rest by an order of magnitude. 
A close look at the eigenvectors corresponding to these 
eigenvalues reveals that these non-affine distortions tend 
to increase the distance between nearest neighbour par¬ 
ticles and reduce next nearest neighbour bond lengths. 
If a nearest-neighbour bond is actually replaced by a 
next-nearest-neighbour one, then the coordination num¬ 
ber of the particles changes and a pair of particles with 
5 and 7 neighbours each would emerge out of the ref¬ 
erence 6-coordinated triangular structure. Each pair of 
neighbouring 5- and 7-coordinated atoms contains a dis¬ 
location (or an anti-dislocation depending on the orien¬ 
tation). These dislocation-anti-dislocation pairs can then 
separate from each other by subsequent non-affine ffuc- 
tuations that change the coordination number of neigh¬ 
bouring atoms. 

Of course in a harmonic lattice defects do not nucleate, 
though non-affine precursor fluctuations exist. Indeed, 
the overlap of particle displacements with a non-affine 
eigenvector b^, given by = b^^A , is a Gaussian 
random variable with probability distribution. 




1 

— , — exp 

y27rb^TCb^ 



an expression analogous to the one for strains. The quan¬ 
tity b^^Cb^ appears as a susceptibility for defect precur¬ 
sor fluctuations. The fluctuation-response relation con¬ 
nects this susceptibility with a response function measur¬ 
ing the response of to a conjugate field. We investi¬ 
gate this connection below. For the rest of this paper, we 
present results only for /i = 1. The corresponding results 
for the degenerate /i = 2 eigenvector are either identical 
or completely analogous. 

To proceed further, we consider the global non-affinity 
X = N~^ Ylf=i averaged over all particles i = 

1... A and introduce a field conjugate to this quantity 
in (§ to obtain, 

H = Hharm-hxNX. (7) 

The extra term in though still quadratic in the par¬ 
ticle coordinates, introduces a many-body force (see sec- 
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FIG. 4. Plot of the inverses of the eigenvalues of the non- 
affine projection PCP together with the eigenvectors for the 
three largest eigenvalues that indicate the most prevalent fluc¬ 
tuations. Note that there are two degenerate soft modes 1,2 
separated by a large gap from the next most important con¬ 
tributor 3. The two soft modes tend to increase the nearest 
neighbour bond length a — b, at the same time decreasing the 
distance between next nearest neighbours c — d (shown by 
gray dashed line). 


tion [V|) that depends on the positions of all particles in 
a given neighbourhood Q. A change in the coordinate 
of particle i modifies not only the local x at that parti¬ 
cle but also those of its neighbours. The force also de¬ 
pends on the reference lattice positions {R}, which act 
as constant parameters. A purely affine transformation 
of O, such as a volume rescaling for example, does not 
produce a non-affine force. This force therefore tracks 
only non-affine distortions away from the reference con¬ 
figuration. The dynamical matrix corresponding to 0 
can be computed without difficulty and therefore the 
statistics of the local x and the local strains, together 
with their space-time correlation functions can be ob¬ 
tained for arbitrary hx us ing the procedures outlined 
in m and sections and This holds true as long 
as the structure of the solid is maintained, i.e. as long as 
the reference configurations {R} remain the global mini¬ 
mum of the modified Hamiltonian 0 . Note that, for our 
MD simulations, in addition to the term proportional to 
hx^ we have also included a small hard core repulsion of 
the Weeks,Chandler, Anderson (WCA) form [16]. This 
prevents atoms from overlapping at large values of hx 
and also introduces anharmonicity in a controlled fash¬ 
ion such that the relative contribution of this term to the 
energy provides us with a measure of anharmonic contri¬ 
butions. The hard core diameter do = 0.6/ was chosen 
to be small enough so that anharmonic effects vanish for 
small values of hx and all our results based on harmonic 
analysis hold in this limit. 

The statistics of the global X and the local x are clearly 
related to each other. For example, the thermal averages 



N-^ 


FIG. 5. Plot of ((AA)^) as a function of N~^ at /ix = 0 from 
our MD simulations for A = 10 x 10,12 x 12,14 x 14,16 x 
16,18 X 18,100 X 100 and 500 x 500 lattices. The straight 
line is the prediction from ^ without any fitting parameters. 
Note that for smaller lattices there are significant deviations 
from the asymptotic slope. 


are equal, {X) = (x). The variance of X is 
((AV)2) = (X2) - {Xf 

R R' 

= iv-'EK x(0)x(R)) - (x)"] 

R 

= Ar-i{(Axf)Eex(R,0) (8) 

R 

The cross-correlation (Ax(R)AA) has the same expres¬ 
sion for any R. 

The variance of X vanishes in the N ^ oo limit as 
expected for an intensive thermodynamic variable (see 
Fig. §; the distribution, P(A), therefore becomes a 
delta function centered at (x)- To obtain the response 
{X{hx))^ for small hx we first compute ((Ax)^) at 
hx = 0 and then use the linear response relation, 

® = ((Ax)^)EC'x(R,0). (9) 

R 

The non-affine field also changes the statistics of the 
dominant displacement fluctuations s^. Since the global 
non-affinity X is quadratic in particle displacements, the 
field hx cannot break the symmetry of 5^. The prob¬ 
ability distribution P{s^) remains Gaussian but with a 
variance (s^) that depends on hx- Again, a linear re¬ 
sponse calculation gives, 

^ = {Xsl)o - (V)o(4)o 

= 2'£hlCPC^h^. (10) 

R 
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FIG. 6. Plots of (a) {X) and (b) (s^), for ji = 1, as a 
function of hx- The points are MD simulation data from 
an A/" = 200 x 200 lattice. Error bars are smaller than the 
size of the symbols. The solid curves through the data are 
the analytic results obtained from direct computation. The 
dash ed straight lines are the linear response predictions ^ X 
(10) respectively. 


with C as given in (|^. 

In Fig. we have plotted {X) and (s^) as functions 
of hx‘ For hx > 0, both of these quantities increase 
making defects more likely to form. In contrast a negative 
hx suppresses those fluctuations that give rise to defects. 
The points are simulation results that are compared with 
the linear response results as well as the full nonlinear 
calculation obtained by evaluating the dynamical matrix 
for the Hamiltonian 0 . 

Space-time correlations of the defect precursors may 
be computed quite straightforwardly from the formalism 
presented in section m Indeed, the correlation of the 
dominant non-affine displacement s^ is (s^)(7s(R, t) = 
(s^(0,0)s^(R,t)) = b^'^Cb^. Fig. shows plots of 
(7^(0, t) and (7s(0,t) against time t for a few values of 
hx- The displacement correlations, like those shown in 
section [TV| are oscillatory and decay slowly in time due 
to destructive interference of the large number of mutu¬ 
ally incommensurate phonon modes that make up these 
localised fluctuations. More importantly. Fig. shows 
that the lifetime of these defect precursors grows as hx 
increases; the time period of the correlation function os¬ 
cillations also increases as expected. Finally in Fig. we 
plot the full Cs(R, t) for three values of hx = 0.00,0.03 
and 0.05 as well as for three values of the time t = 0, 2 and 
5 as in Figs, [^an di Unlike the correlation functions for 
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FIG. 7. Plots of the normalised (a) Cx(0,t) and (b) Cs(0,t) 
for three values of hx compared with results from MD simu¬ 
lations of 500 X 500 particles. Note that the precursor fluctu¬ 
ations become more long-lived as hx increases. 


X and 5 , the correlations of are “anti-ferromagnetic”, 
i.e. a fluctuation 5 ^ of any sign at some lattice point 
induces a fluctuation of s^ of the opposite sign at the 
neighbouring lattice point. 

What is the effect of external stress 5] on the defect 
precursors? It is again easy to answer this question by 
a straightforward calculation: one only needs to include 
the term ^Fe Hamiltonian ([^, where 

e(R i) is now the loeal strain at particle i. The probability 
P(s^) remains Gaussian with the same variance but now 
the ± symmetry of is explicitly broken and P(s^) 
is shifted with a mean {s^) = 7 ^ 0. To 

lowest order, therefore - or exactly in a harmonic solid - 
stress biases the distribution of defect precursors without 
changing their variance. Similarly, I] does not affect the 
space-time correlation functions of s^. 


V. GENERATING hx USING LASER 
TWEEZERS 


In this section, we propose an experimental realisa¬ 
tion of the many-body term in the Hamiltonian 0 using 
dynamic laser traps - a technology currently available 
within most sophisticated experimental optics research 
groups. Colloidal particles are dielectric and therefore 
become polarised in an electric field. Fairly intense light 
from a laser may be used to trap these particles, which 
experience a force proportional to the gradient of the 
light intensity /(r) and therefore prefer to accumulate 
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hx — 0.03 



FIG. 8. Plot of the normalised Cs(R, t) for hx = 0.00,0.03 and 0.05 for t = 0, 2 and 5. Note that the precursor fluctuations 
are spatially anisotropic and anti-correlated at all times. The non-affine held makes the correlations stronger. We have plotted 
the correlations for /x = 1; the corresponding functions for // = 2 are essentially rotated by 90°. 


in regions of large /(r). This effect is extremely useful 
in manipulating colloidal beads in the lab to investigate 
myriads of phenomena from biology to material science. 
There are many reviews and books on the subject, such 
as Ref. [13]. More specihcally, optical traps have been 
used to manipulate colloidal solids, introduce defects and 
watch their dynamics using video microscopy m- 

The term proportional to hx in 0 depends only on 
the reference lattice set {R} and the instantaneous par¬ 
ticle positions and can be generated for every particle i 
once a particular conhguration is known. For example, 
one can explicitly write, 

Xi = - Ui)^Pj-i,k-i{uk - Uj) 

jk 

where we gather the cartesian components of P for a given 
pair of particles into a matrix and assume that this 

matrix is zero when j or k are outside the neighbourhood 
Q around i. Then 

NX = - u*) 

ijk 

The force on particle i is = —{d/dri){—hxNX). Di¬ 


rect differentiation of the expression for X then gives 
F^ = ‘2ihx ^ ^ [Fj—i,/c—iu/c) T Pi—ji (u/c w)] 

jk 

( 11 ) 

The first contribution comes from Xii fhe second from Xj 
with j i. 

The fact that the above forces can be worked out from 
the positions of the particles and their nearby neigh¬ 
bours (nearest and next-nearest neighbours, if the coarse- 
graining volume Q contains exactly the nearest neigh¬ 
bours) suggests the following algorithm for generating a 
uniform non-afiine field hx for a set of N colloidal par¬ 
ticles: 

1. At any instant obtain the coordinates of the N par¬ 
ticles through video microscopy. 

2. Randomly choose a subset of M of these particles 
that will have /ix-forces applied to them. 

3. For each of the M particles obtain the values of 
the necessary forces from the coordinates of their 
neighbors. 
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4. Apply the forces by constructing a set of M laser 
traps. The traps will need to be placed slightly 
away from the respective present particle positions 
so that the particles experience exactly the forces 
calculated from (11). The exact displacements of 


the traps will depend on /(r) and therefore vary 
with the specific apparatus and implementation. 

5. In the next instant repeat steps 1 — 4 above, choos¬ 
ing another random subset of M particles to track. 


If these steps are repeated on a time scale much faster 
than the typical diffusion time of colloids, then one should 
be able to simulate a uniform field hx applied across all 
the N particles. It is possible to update dynamical traps 
at 200 — 600 Hz, and set up at least M = 300 traps simul¬ 
taneously for micron sized colloidal particles using spatial 
light modulator (SLM) technology [13]. This should be 
enough to generate a uniform hx as long as the ratio of 
the dynamical timescale to the update timescale is larger 
than N/M. Alternatively, one may also look at the effect 
of a local hx which couples to the x of a single particle 
and can create local defect precursors. Statistics of such 
local and non-uniform, dynamic, light fields may also be 
computed, if desired, from the formalism outlined in this 
work. 


VI. DISCUSSION AND CONCLUSIONS 

In this paper we have calculated the space-time correla¬ 
tion functions for thermally generated non-afhne fluctua¬ 
tions and elastic strains in a harmonic ideal crystal. The 
non-affine and elastic strain fields were obtained by pro¬ 
jecting atomic displacements into orthogonal afhne and 
non-affine sub-spaces defined by coarse-graining over a 
fixed volume U. Our results show that these correlation 
functions decay to zero with time and over distance al¬ 
though the relaxation to the late time value is oscillatory 
rather than monotonic. The time correlation functions 
for non-afhne huctuations and strains have not been de¬ 
scribed so far in the literature though we feel that they 
may be obtained easily for colloidal solids using video 
microscopy. This should allow verihcation of our re¬ 
sults against experimental data BED]. Note that the 
harmonic approximation that we have used throughout 
has been demonstrated to describe colloidal solids rather 
well [9]. 

In addition we have identihed particular non-affine 
huctuations in the 2 d triangular lattice which, we demon¬ 
strate, are precursors to the production of dislocation- 
anti-dislocation pairs and arise naturally from a system¬ 
atic coarse-graining procedure. We emphasise that the 
defect precursors 5 ^ are not themselves defects since the 
equilibrium average (s^) = 0 . 

In order to form dislocation pairs, these localised huc¬ 
tuations need to condense by escaping over a, possibly 
stress dependent, barrier A/, a process not describable 
within harmonic theory [35]. Indeed, if the bond c-d in 


Fig. a were to form, a Burgers circuit around particle 0 
would yield a non-zero Burgers vector. One can argue, as 
below, that the non-affine held hx will actually greatly 
enhance the formation of such dislocation dipoles in a 
real solid. The rate of barrier crossing is proportional 
to exp(—;dA/) with a prefactor, the so called “attempt 
frequency” which is a product of the characteristic fre¬ 
quencies of oscillation of the system in its parent state 
and at the saddle point [35]. Consider the neighbour¬ 
hood U of a single particle. In a solid with anharmonic 
forces between particles, the free energy for producing 
a precursor huctuation of amplitude has the form 
f{s^) = —where A, B and C are, possibly 

temperature (and stress) dependent, phenomenological 
parameters. Note that, in the harmonic limit H = C = 0 
and A oc ( 5 ^)“^. This form for the free energy ensures 
that the symmetry is preserved and a non-zero bar¬ 
rier for the nucleation of a dislocation dipole ((s^) 7 ^ 0 ), 
given by the saddle point value of /(s^), exists. When 
hx is turned on, this has the effect of increasing (s^) 
(see ([To])). This has two consequences: it decreases both 
the attempt frequency and A/ with the latter effect far 
outweighing the former and effectively causing an overall 
increase in the rate of production of dislocation dipoles. 
For negative hx, on the other hand A/ is increased and 
dislocation nucleation is suppressed. 

The dynamics considered in our formulation is entirely 
composed of lattice vibrations. In a crystalline solid one 
needs to consider, in addition, the slow vacancy diffu¬ 
sion mode [36] . Since the vacancy concentration in crys¬ 
talline solids at temperatures far from melting is vanish¬ 
ingly small, this contribution is mostly negligible at low 
temperatures. However, close to the melting transition, 
the diffusion of vacancies does contribute significantly. 
Within a harmonic theory, there is, of course, no descrip¬ 
tion of vacancy diffusion. On the other hand, vacancy 
diffusion over large distances occurs by small movements 
of atoms across distances of the order of the lattice spac¬ 
ing — and so of the sort involved in the nucleation of 
a dislocation pair. Hence precursor fluctuations for va¬ 
cancy diffusion may be similarly described as a non-affine 
distortion of a volume U containing a single vacancy. A 
calculation of vacancy migration precursors using a pro¬ 
cedure similar to the one described in this work is in 
progress. 

Our calculations may also be generalised to amorphous 
solids. In such solids, the lack of a clearly defined ref¬ 
erence configuration makes the identity of the relevant 
non-affine fluctuations debatable. The dominant defor¬ 
mation mechanisms in amorphous solids are local atomic 
rearrangements that resemble, somewhat, our defect pre¬ 
cursor modes. However, it is not clear whether, in a 
particular realisation of the amorphous structure, such 
precursors are “frozen in” laDlEl] even at zero stress or 
are produced during the deformation protocol [32]. We 
believe that a generalisation of our calculation may be 
able to elucidate this point by looking at neighbourhoods 
with large non-affine susceptibility and determining the 
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response to both stress and the non-affine field hx- 
The results presented here may be verified in experi¬ 
ments on colloidal solids or dusty plasmas in the presence 
of a field hx produced using laser tweezers [iHiin]. Since 
functionalised colloidal assemblies have many technolog¬ 
ical applications, control over their structure may be of 
some use [12]. In section |V| we outline an algorithm 
which, we believe, can be implemented in practice. Sim¬ 
ilar ideas have been reported in the literature where 
light fields have been used to create dislocations and grain 
boundaries by manipulating individual colloidal parti¬ 
cles. We believe our approach allows greater control by 
targeting, instead, defect precursor fluctuations. First 
of all, one is able to both increase as well as suppress 
defect densities in a crystal by an external light held. 
Also, if hx is applied sufficiently slowly, the solid may be 
persuaded to remain in thermodynamic equilibrium at a 
given temperature throughout the process without pro¬ 
ducing unwanted stresses and deformations. Finally, the 
specific dynamics of such protocols (switching hx off or 
on at some rate) can be computed within the formalism 
discussed here. It is also, in principle, possible to excite a 
local non-afflne displacement or even a specific non-afflne 
mode, say si, at a specific point using our ideas. For the 
latter case, however, one needs to know beforehand the 
eigenvectors of the local PCP, which involves a knowl¬ 


edge of the interactions embodied in the C matrix. This 
introduces uncertainties that are not encountered while 
imposing hx- For dusty plasmas [12], the equations we 
have used for the space-time correlation functions are im¬ 
mediately applicable. For colloidal particles dispersed in 
a liquid, of course, one needs to account for damping 
and Brownian noise terms in the dynamical equation @ 
to compare time-dependent quantities with experiments. 
Equilibrium predictions, though, would continue to be 
valid. Also anharmonic interactions, always present in 
real colloids, would lead to metastable defects at positive 
hx- Our calculations are then directly valid for small 
values of the held before such nucleation events actually 
take place. We believe that in this case, our results will 
be of much value for checking and validating the relevant 
experiments. 
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